Physica A 354 (2005) 249-261 



Fractional Ginzburg-Landau equation 
for fractal media 

Vasily E. Tarasov 1 and George M. Zaslavsky 2 ' 3 

1) Skobeltsyn Institute of Nuclear Physics, 
Moscow State University, Moscow 119992, Russia 
2) Courant Institute of Mathematical Sciences, New York University 
251 Mercer St., New York, NY 10012, USA, 
3) Department of Physics, New York University, 
2-4 Washington Place, New York, NY 10003, USA 



Abstract 

We derive the fractional generalization of the Ginzburg-Landau equation from the 
variational Euler-Lagrange equation for fractal media. To describe fractal media we 
use the fractional integrals considered as approximations of integrals on fractals. Some 
simple solutions of the Ginzburg-Landau equation for fractal media are considered and 
different forms of the fractional Ginzburg-Landau equation or nonlinear Schrodinger 
equation with fractional derivatives are presented. The Agrawal variational principle 
and its generalization have been applied. 
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1 Introduction 



We call the fractional equation an equation that contains fractional derivatives or integrals. 
Derivatives and integrals of fractional order have found many applications in recent studies 
in physics. The interest in fractional equations has been growing continually during the 
last few years because of numerous applications: kinetic theories of systems with chaotic 
dynamics [UI2]; pseudochaotic dynamics [3]; dynamics in a complex or porous media [3]; 
random walks with a memory and flights El El ; and many others. 

Fractional integrals and derivatives can be used to describe processes with different levels 
of complexity. The new type of problems has rapidly increased interest in areas in which 
the fractal features of a process or the medium impose a necessity of applying tools that are 
non-traditional in "regular" smooth physical equations. In many problems the real fractal 
structure of matter can be disregarded and the medium can be replaced by some smoothed 
continuous mathematical model where fractional integrals appear ||. The order of the 
fractional integral is equal to the fractal mass dimension of the medium and in this way one 
can take into account the fractality of the media. It was proved that integrals on a net of 
fractals can be approximated by fractional integrals [5] and that fractional integrals can be 
considered as integrals over the space with fractional dimension up to a numerical factor 
[TU] . The last proof used the formulas of dimensional regularizations [TTj . 

The fractional generalization of the Ginzburg-Landau equation was suggested in Ref. 
[T2] . This equation can be used to describe the dynamical processes in a media with fractal 
dispersion. Since the fractals can be realized in nature as a fractal process or fractal media, 
it is interesting to derive fractional Ginzburg-Landau (FGL) equation using a corresponding 
generalization of the free energy functional. 

It is known ^3] that the Ginzburg-Landau equation 



gAZ - aZ - bZ 3 = 



can be derived as the variational Euler-Lagrange equation 



8F{Z(x)} 
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for the free energy functional 

F{Z(x)} = F + i / [^(VZ) 2 + aZ 2 + b -Z 4 }dV 3 , (2) 

where the integration is over three-dimensional region W of continuous media. Here F is a 
free energy of the normal state, i.e. F{Z} for Z = 0. In this paper we consider a fractional 
generalization of (J2J) that can appear from two places: fractional generalization of the integral 
in Eq. (J2J) and fractional generalization of the derivatives in Eq. (J2J). Different forms of the 
FGL equation are considered in the coordinate space and in the dual space after applying 
the Fourier transform. Finally, the Agrawal variational principle [2] and its generalization 
have been applied to obtain the FGL equation. 



2 Fractional generalization of free energy functional 

Let us consider the thermodynamic potential (free energy functional) F{Z(x)} for the non- 
equilibrium state of the fractal medium, where Z(x) is a scalar field. The free energy 
functional has the form 

F{Z(x)} = F + / .F(Z(x), VZ(x))dV 3 , (3) 
Jw 

where JF(Z(x), VZ(x)) is the free energy density; VZ = dZ/dx. For the Ginzburg-Landau 
potential (J2J this density is 

^(Z(x), VZ(x)) = \[g{VZf + aZ 2 + ^Z 4 ]. 

The simplest fractional generalization of (J2J) can be written in the form 

F{Z(x)} = F + [ F{Z(x), VZ(*))dV D . (4) 

Here D is a fractal mass dimension of the fractal medium, and dVu is an element of the 
D-dimensional volume: 

dV D = C 3 (D,x)dV 3 . (5) 

Note that the interpretation of the fractional integration is connected with fractional dimen- 
sion PU], which follows from the well-known formulas for dimensional regularizations [TT] : 



the fractional integral can be considered as an integral in the fractional dimension space up 
to the numerical factor T(D/2)/(2tt d ^T(D)). 

For the Riesz definition of fractional integral, the function Cs(D, x) is 

3V ' ; r(D/2) 11 v ; 

The initial point of the fractional integral is set to zero and the numerical factor in Eq. 
(JHJ) gives the usual integral in the limit D — > (3 — 0). Note that the numerical factor 
%\D) = r(l/2)/2 D 7r 3 / 2 r( J D/2) used in Ref. £] leads to 7^(3-0) = r(l/2)/2 3 7r 3 / 2 r(3/2) 
in the limit D -»• (3 - 0). 

For the Riemann-Liouville fractional integral, the function Cs(D,x) is 

ft(Ax) = J5^^. (7) 

In Cartesian's coordinates x^. 

3 

x = ^Ja^e^, = dx\dx2<ix^. 

k=l 

The simplest fractional generalization of (J2J) for fractal media is 

F{Z(x)} = F + l [ [g(VZy + aZ 2 + V]dVb. (8) 

Equation (jHJ) can be considered as terms of the expansion of a functional in series over small 
values Z(x) and the integer derivatives VZ = dZ/dx.. 

For homogeneous media without external fields the parameter Z does not depend on 
coordinates and Eq. (jHJ) gives 

F{Z} = F + ^ + b -f Z \ (9) 

where Vd is a D-dimensional volume of fractal medium in the region W. The equilibrium 
value of Z corresponds to the minimum condition for 

dF{Z}/dZ = 0. 

If a/6 > 0, then the free energy potential has the single minimum Z = 0. If a/6 < 0, 



then there are two minima Z = ±y^—a/b. The phase transition is realized when a/6 passes 
through zero. 



3 Ginzburg-Landau equation for fractal media 

In the general case, the equilibrium value of Z(x) is defined by the condition of the minimum 
of F{Z(x)}, which has the form of variational Euler-Lagrange equation (0). Using functional 
(jlj with the free energy density J(Z(x),VZ(x)), we get the Euler-Lagrange equation 




where C 3 (D, x) is defined by Eq. (JHJ) or (JJJ), and = d/dx k , k = 1,2, 3. Using (JSJ), we get 
the equation 

gC^(D,^)V k (c 3 (D,x.)V k Z^ -aZ-bZ 3 = (11) 

that can be considered as the Ginzburg-Landau equation for fractal media in the case T = 
^(Z(x),VZ(x)). 

Eq. (|TT|) can be rewritten in an equivalent form 

gAZ + E k (D,x)V k Z - aZ -bZ 3 = 0, (12) 

where 

E k (D,x) = C 3 - 1 ( J D,x)V fc C 3 p,x). 
In one-dimensional case, when < D < 1, we have 

^ + - Ca/P)^ - (W^ 3 = 0, (13) 

where x = x\, and Z' x = dZ/dx\. The simplest solutions of this equation are considered in 
Section 7. Note that this equation has an analog in the form of an equation for a nonlinear 
oscillator with friction. Therefore, the Ginzburg-Landau equation for fractal media describes 
dissipative nonlinear oscillations in coordinate space. 

In the general case, the free energy density can also depend on fractional derivatives of 
Z. In this case we can use the fractional generalization of the Agrawal variational equation 
[T3] . This case will be considered in Section 6. 
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4 Ginzburg-Landau equation with fractional derivatives 

Let us recall the appearance of the Ginzburg-Landau equation with fractional derivatives 
suggested in Ref. [12]. Consider wave propagation in some media and present the wave 
vector k in the form 

k = k + K = k + K|| + (14) 

where k is the unperturbed wave vector and subscripts (||, _L) are taken respectively to the 
direction of ko. Considering a symmetric dispersion law uo = u(k) for wave propagation with 
k /eg, we have 

u(k) = co(\k + k\) uj(k ) + c(|k + k\ — k ) « co(k ) + ck\\ + rr-Kj., (15) 

where c = du/dk . Expression (|T3jl in dual space ("momentum representation") corresponds 
to the following equation in coordinate space 

dZ dZ c . „ 

-« = K &7 + ^ AZ (16) 

with respect to the field Z = Z(x, t), where x\ is along k . Here we use the operator 
correspondence between the dual space and the usual space-time: 

*» ~ ( ^ )2 — -^ = ~h~h (17) 

A generalization to the nonlinear case can be carried out as in ([T3]l through a nonlinear 
dispersion law dependence on the wave amplitude: 

to = u(k, \Z\ 2 ) « u(k,0) + b\Z\ 2 = u(\k o + K\,0) + b\Z\ 2 (18) 

with some constant b = du(k, \Z\ 2 )/d\Z\ 2 at \Z\ 2 = 0. In analogy with (|16p. the nonlinear 
equation takes the form 

87 r)7 r 

This equation is also known as the nonlinear Schrodinger (NLS) equation in which any of 
the coefficients may be complex. Indeed, for Z = Z(t,X\ — ct, x 2 ,x 3 ) we have 

r)7 r 

-i-K7 = ^r^z-^z-b\z\ 2 z (20) 

at 2k 



with u = u(k ). 

Wave propagation in media with fractal properties can be easily generalized by rewriting 
the dispersion law (|15j). Namely, one can replace (JT5j) and (fTHj) by the following equation in 
dual space: 

u{k, \Z\ 2 ) = cu(k , 0) + ck\\ + c a (K 2 ± ) a/2 + b\Z\ 2 (21) 

with a fractional value of 1 < a < 2 and a new constant c a . 

Using the connection between fractional derivative and its Fourier transform ^5] 

(_A)«/ 2 < — > (Kl) a/2 (22) 

we obtain the equation 

=i%-%.-W»Z-<*Z-Wz, (23) 

which corresponds to (|21j) in time-coordinate space. For Z(t,x% — ct,X2,Xs) traveling waves 
as in (j20J), we get the equation 

^ = ^(-&T /2 Z + oo Z + b\Z\ 2 Z (24) 
at 2k 

that can be called nonstationary FGL equation or fractional nonlinear Schrodinger equation 
(FNLS). Let us comment on the physical structure of (|24|) . The first term on the right-hand 
side is related to the wave propagation in a media with fractal properties. The fractional 
derivative may also appear as a result of ray chaos [THJ HZj or due to super diffusive wave 
propagation (see also the discussion in Refs. [TJ HHj and corresponding references therein). 
Other terms on the right-hand-side of Eqs. (|2Hjl . and ()24j) correspond to wave interaction 
due to the nonlinear properties of the media. Thus, Eq. (J2H) can describe fractal processes 
of self-focusing and related issues. 

We may consider a one- dimensional simplification of (J24|) . i.e. 

dZ 

c— = gB*Z + aZ + b\Z\ 2 Z (25) 
with some generalized constants g,b,c, and the Riesz fractional derivative 

(d«/x*)= - — (f mdz + rtiim^) (26) 



where n — 1 < a < n. As a particular case, we can reduce Eq. (f23j) to a propagating wave 
solution Z = Z(x), x = x% — ct. In this case, equation ()25|) takes the form of the stationary 
FGL: 

gB°Z + cDlZ + aZ + b\Z\ 2 Z = (27) 

or for real Z: 

gV*Z + cQ\Z + aZ + bZ 3 = 0. (28) 



Eq. (j28|) differs from the fractional Burgers equation ^81 EI] in the structure of the nonlinear 
term. Nevertheless, an analysis similar to that in Refs. ^Hl EH] may be performed to obtain 
some estimates on the solution. 

It is well known that the nonlinear term in equations of the type (I2UJ) leads to a steepening 
of the solution and its singularity. The steepening process may be stopped by a diffusional 
or dispersional term, i.e. by a higher derivative term. A similar phenomenon may appear 
for the fractional nonlinear equations (|23j) . (}2Tj) . It has been shown in Ref. [T^] that for the 
fractional Burgers equation there exists a critical value a c such that the solution is regular 
for all x if a > a c . 

5 Ginzburg-Landau equations with fractional deriva- 
tives in dual space 

In order to derive the FGL equations from a variational principle, we can use two methods. 
One of them was suggested in Ref. ^U] and it will be considered here. The second one is 
a generalization of the Agrawal variational equations for the fractal media, and will be 
considered in Section 6. 

Let us introduce the derivatives with a fractional power of coordinates : 

d Ixkl 1 '^ d 



dxl P dx k 



(29) 



These derivatives transform into fractional derivatives in dual space ("momentum represen- 
tation"). The Fourier transform of the coordinates |xfc| 1_/3 leads us to the Reisz fractional 



derivative 



\Xk 



1-/3 _ D l-/9 



Ilk 



where we use definition (|26|) . The Fourier transform of Eq. (|29|) is 



9 W"" 9 I D 'A, (30) 



Consider free energy functional T = T{Z, dZ/dx^ ) for one-dimensional case. Then the 
FGL equation for fractal media is of the form 

gCt\P,x)^(c 1 tf,x)-^z)-aZ-bZ* = U, (31) 

where Ci((3,x) = or in an equivalent form 

A N ^^ Z _ aZ _ 6Z » = . (32 ) 

After Fourier transform of this equation, we obtain the FGL equation in dual space: 

^p\- 2 \y 2 Z) + aZ + bZ 3 = 0. (33) 

Here Z(y) is a Fourier transform of Z(x), i.e. Z(x) < — ► Z(y). For the field ^f(y) = y 2 Z(y), 
we can rewrite Eq. ()33)1 in the form 

gT) a x m(x) + a(x)V(x) + b{x)^ 3 {x) = 0, 

where a = 2 - 2/3, < a < 2 and a(x) = a(3 2 /\x\ 2 , b(x) = b/3 2 /\x\ 6 . 

As a result, the FGL equation for fractal media has fractional derivatives in dual space. 
The FGL equation with fractional derivatives in coordinate space is considered in the next 
section. 

6 Ginzburg-Landau equations with fractional deriva- 
tives from fractional Agrawal principle 



Let us, first, extend the fractional Agrawal variational equations [M] for the case of fractal 
media. 
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The simplest generalization of the free energy functional with fractional integrals and 
derivatives can be written in the form 



Jw 



(34) 



Here D is a fractal mass dimension of fractal medium, and dVn is an element of D- 
dimensional volume such that dV D = C 3 (D,^.)dV 3 . The function C 3 (.D,x) is defined by 
relations © and (J7|). D a is a fractional derivative that is defined by Eq. (J2HJ). The potential 
density is 

F= 1 -[g(B a Zr + aZ*+ h -Z% (35) 

The condition of the minimum of free energy potential ()34|). (|33j) gives the Euler-Lagrange 
equation in the form 

C 3(Ax)g + X>: fc (c 3 (Ax)^) =0. (36) 

fc = l X k 

In general, D ^ 3a and Eq. (j3T?j) gives 

3 

gC^(D, x) £ D« t (c 3 (A x)D« fc Z) + + 6Z 3 = 0. (37) 
fc=i 

It can be considered as a fractional generalization of the Ginzburg-Landau equation for 
fractal media. Consider some cases of this equation. 

1. In the one-dimensional case Z = Z(x), the coordinate fractional derivative is D", i.e., 

T = F{Z,T> a x Z). 

The FGL potential density has the form 



t =\ g {p-z)\\z^\z\ (at 



Using the formulas for fractional integration by parts 

ft r < 

/oo poo 
f(x)B a x g(x)dx = / g(x)B a J(x)dx, (39) 
-oo J — oo 
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we obtain the Euler-Lagrange equation 

D " i Cl ( A x) d5£z) + Cl ( A x) % = °' (40) 

where the function C\(D,x) is defined as 



\x 



D-l 



C 1 (D,x) = ^—-. (41) 
Using for T (jHHj) . we arrive at 

Cf 1 ( J D,x)D^C 1 (,D,x)D^) +aZ + 6Z 3 = 0. (42) 
For the case D = 1, we have C*i = 1 and (|41j) transforms into 

D 2a Z + aZ + 6Z 3 = 0, (43) 

where is the Riesz derivative. 

2. The free energy functional density T can be a function of Z = Z(x), and two different 
fractional derivatives D"Z, and DfZ, a ^ /3, i.e., 

^ = ^(Z,D^,DfZ). 

The FGL potential density has the form 

^ = ^i(D^) 2 + i, 2 (DfZ) 2 + \z* + \z\ (44) 

For example, a = l,0</3<l. Using again (J59j) . we obtain the following Euler-Lagrange 
equation, 

< c ^^z) + < c ^^) + = °- (45) 

where Ci(D,x) is given in (|41jl. For the FGL density (|44j). we get the FGL equation for 
fractal media 

^iCf 1 ( J D,x)D^Ci( J D,x)D^) +^ 2 Cr 1 ( J D,x)Df(c 1 (D,x)Dfz) + aZ + bZ 3 = 0. (46) 
Particularly, when D = 1 and G\ = 1, we have 

^iDfZ + £ 2 D 2/3 Z + aZ + 6Z 3 = 0, < a, < 1. 
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3. In the general case, the free energy density functional is a function of Z = Z(x), and 
derivatives D"*Z of fractional orders different for coordinates Xk, i.e., 



T = F(Z, T>%Z, T>T 2 Z, D°fZ). (47) 
For this case, the FGL equation for fractal media is 

3 

gC^(D, x) D S (C 3 (D, x)Dgz) + aZ + bZ 3 = 0. (48) 
fc=i 

The sum of orders can be equal to the fractal mass dimension D of the medium 

cki + a 2 + a 3 = D, 
but in the general case it can be that «i + «2 + «3 7^ D. 



7 Simplest solutions 



Let us consider a simple example of the Ginzburg-Landau equation ([lip for fractal media in 
the one-dimentional case. Then 

dZ dZ 

and Z = Z(xi). Eq. (|31j) for one- dimensional case f}\ = D and < D < 1 transforms into 

gC^{D,x) — (c x {D,x)—Z) - aZ - bZ 3 = 0, (49) 

where x = Xi, and Ci(D,x) is defined by (jJTJ). This equation can be rewritten in an 
equivalent form 

Z'L + ^Z' x - (a/g)Z - (b/g)Z* = 0, (50) 

9m 

where Z' x = dZ(x)/dx. Note that Eq. ([5()p has an analog of the equation for a nonlinear 
oscillator with friction 

x{t) = ax(t) + (3x 3 (t) - iit)x{t), 

where we use 

a = a/g, (3 = b/g, 7 (t) = (D - l)/g\t\. 
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Therefore the FGL equation for fractal media describes coordinate nonlinear oscillations 
with a "dissipative"-like term. 

Let us consider solution of with 6 = 0. Z(x) satisfies the equation 

#[Ci(A x)Z'J x - ad(D, x)Z = 0, 

where x G (0; oo), that can be rewritten as 

gxZ'l x {x) + (D - l)Z' x (x) - axZ(x) = 0. (51) 

The corresponding solution is 

Z(x) = C l x l - D l 2 J v {^aJgx)+C 2 x 1 - D/2 Y v {^aJgx), 

where v = |1 — D/2\, J v {x), and Y v (x) are the Bessel functions of the first and second kind. 
The solution of Eq. (f5T|) with the conditions 

Z(l) = e, Z'(l) = c 

has the form 

Z(x) = a; 1 - D/2 J D /2-i(v /Z Ax)(c J Dv /Z A Y D/2 {V^A) - {eA)^A Y D/2 (V=A)+ 
+cAY D/2+1 (V^A))/(AV^A( Jd/2+i(V=A) Y d/2 (V^A) - Y d/2+1 (V^A) J d/2 (v^4))) + 

+x 1 - D/2 Y D/2 „ 1 (V^Ax)(-cD J D/2 (V^A) + ^Ac J D/2+1 (V^A)+ 
+eA J d/2 (V^A))/(a( J D/2+1 (V^A) Y d/2 (V^A) - Y D/2+1 (V^A) J D/2 (V=A)fj 
where A = a/g. 

The solutions with b — 0, g — 1, and a/g = — 1 are demonstrated in Fig. 1 for D = 9/10 
and in Fig. 2 for D = 1/10. The growth of solution is saturated by the nonlinear term and 
full solution for (|50|) may appear similar to the vicinity of a stable limit cycle. 
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Figure 1: Solution for b = 0, a/g = — 1 and D = 9/10. 
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8 Conclusion 



Derivatives and integrals of fractional order have found many applications in recent stud- 
ies in physics. In this paper, we use the fractional integrals to describe fractal media with 
noninteger mass dimensions. Fractional integration can also be used to describe time depen- 
dent processes in fractal media. The fractional Ginzburg-Landau equation for fractal media 
is derived from the corresponding generalization of free energy functional and variational 
Euler-Lagrange equations. Generalization of the Agrawal variational equation 14 a appears 
to be very useful for this goal. 

The fractional continuous models have wide applications. Partly, it is because relatively 
small number of parameters can define a complex medium with a rich structure. In many 
cases, the real fractal structure of matter can be disregarded and the medium can be de- 
scribed by some fractional continuous mathematical model. Smoothing of the microscopic 
characteristics over the physically infinitesimal volume transform the initial fractal medium 
into a fractional continuous model that uses the fractional integrals of the same order as the 
fractal mass dimension of the media. The fractional integrals can be considered as an ap- 
proximation of the integrals on fractals |9| . Note that fractional integrals can be considered 
as integrals over the space with fractional dimension up to a numerical factor ^U|. Some 
applications of the FGL equation was discussed also in Ref. [2*U| . 
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